QCD corrections to polarization of J/ip and T at Tevatron and LHC 
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In this work, we present more details of the calculation on the next to leading order (NLO) 
QCD corrections to polarization of direct J/ip production via color singlet at Tevatron and LHC, 
together with the results for T for the first time. Our results show that the J/^ polarization status 
drastically changes from transverse polarization dominant at leading order (LO) into longitudinal 
polarization dominant in the whole range of the transverse momentum pt of J/ip when the NLO 
corrections are counted. For T production, the pt distribution of the polarization status behaves 
almost the same as that for J/ip except that the NLO result is transverse polarization at small 
Pt range. Although the theoretical evaluation predicts a larger longitudinal polarization than the 
measured value at Tevatron, it may provide a solution towards the previous large discrepancy for 
J/ip and T polarization between theoretical prediction and experimental measurement, and suggests 
that the next important step is to calculate the NLO corrections to hadronproduction of color octet 
state J/-!/)*-*' and T'-^^ Our calculations are performed in two ways, namely we do and do not 
analytically sum over the polarizations, and then check them with each other. 

PACS numbers: 12.38.Bx, 13.25.Gv, 13.60.Lc 



I. INTRODUCTION 



The study of J/ip production on various experiments is 
a very interesting topic since its discovery in 1974. It is a 
good place to probe both perturbative and nonperturba- 
tive aspects of QCD dynamics. To describe the huge dis- 
crepancy of the high-pt J/ ip production between the the- 
oretical calculation based on color singlet mechanism [l| 
and the experimental measurement by the CDF collab- 
oration at the Tevatron^, color-octet mechanismQ was 
proposed based on the non-relativistic QCD(NRQCD)[|. 
The factorization formalism of NRQCD provides a theo- 
retical framework to the treatment of heavy-quarkonium 
production. It allows consistent theoretical prediction to 
be made and to be improved systematically in the QCD 
coupling constant and the heavy-quark relative ve- 
locity V. The color singlet mechanism is straightforward 
from the perturbative QCD, but the color-octet mecha- 
nism depends on nonperturbative universal NRQCD ma- 
trix elements. So various efforts have been made to con- 
firm this mechanism, or to fix the magnitudes of the 
universal NRQCD matrix elements. Although it seems 
to show qualitative agreements with experimental data, 
there are certain difficulties in the quantitative estimate 
in NRQCD for J/ih photoproduction at the DESY ep col- 
lider HERA J/^iV) and T polarization 
of hadronproduction at the Fermilab Tevatron, and J/ip 
production in B-factories. A review of the situation could 
be found in Ref. |ll|. 

Without NLO corrections, it is difficult to obtain 
agreement between the experimental results and leading 
order theoretical predictions for J/ip production. There 
are a few examples shown that NLO corrections are quite 
large. It was found that the current experimental re- 
sults on inelastic J/ip photoproduction[l2, E^l are ade- 
quately described by the color singlet channel alone once 



higher-order QCD corrections are included B, [1]. Al- 
though ref. ^14,] found that the DELPHI [ij data evi- 
dently favor the NRQCD formalism for J/ip production 
7-I-7 J/ip + X, rather than the color-singlet model. 
And it was also found in ref. [lB| that the QCD higher or- 
der process 7 + 7^ J/iJj + c + c gives the same order and 
even larger contribution at high pt than the leading or- 
der color singlet processes. In reL [l^, the NLO process 
c + g J/ip-\-c where the initial c quark is the intrinsic c 
quark from proton at Tevatron, gives larger contribution 
at high Pt than the leading order color singlet processes. 
The large discrepancies found in the single and double 
charmonium production in e^e^ annihilation at B facto- 
ries between LO theoretical predictions [H, [l^ [13] and 
experimental results [21I [2^ were studied in many work. 
It seems that they may be resolved by including higher 
order correction: NLO QCD and relativistic corrections 
[ll[2l[2i,[2i,[2i,[27|,i28,]. 



Based on NRQCD, the LO calculation predicts a siz- 
able transverse polarization for J/ip production at high 
Pt at Tevatron [23J33, [Ml while the measurement at Fer- 
milab Tevatron [32| gives slight longitudinal polarized 
result. In a recent paper [33], the measurement on po- 
larization of T production at Tevatron is presented and 
the NRQCD predication [s^] is not coincide with it. Be- 
yond the NRQCD framework, there is a try by using s- 
channel treatment to J/ip hadronproduction in the work 
of ref [3^, which gives longitudinal polarization. Within 
the NRQCD framework, to calculate higher order correc- 
tions is an important step towards the solution of such 
puzzles. Recently, NLO QCD corrections to J/ip hadron- 
production have been calculated in ref [s^. The results 
show that the total cross section is boosted by a factor of 
about 2 and the J/ip transverse momentum pt distribu- 
tion is enhanced more and more as pt becomes larger. A 
real correction process g + g — > J/ip + c+c at NLO, which 
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is not included in the ref. [s^ , was calculated in [T^ [sj ■ 
It gives sizable contribution to pt distribution of J/ ip at 
high Pt region, and it alone gives almost unpolarized re- 
sult. Therefore it is very interesting to know the result 
of J/ip polarization when NLO QCD corrections are in- 
cluded. In a recent Letter [sl] , we presented a calculation 
on the NLO QCD corrections to the J/ip polarization in 
hadronproduction at Tevatron and LHC. In this paper, 
we give more details of the calculation, and the results for 
T polarization for the first time. The results show that 
the polarizations of J/ijj and T are drastically changed 
from more transverse polarization at LO into more longi- 
tudinal polarization at NLO. Meanwhile, our results for 
total cross section and transverse momentum distribution 
is consistent with ref. ^36*1 . In this calculation, we use our 
Feynman Diagram Calculation package (FDC)[3^ with 
newly added part of a complete set of method to calculate 
tensor and scalar integrals in dimensional regularization, 
which was used in our previous work[24, 25]. 

This paper is organized as follows. In Sec. II, we give 
the LO cross section for the process. The calculation of 
NLO QCD corrections are described in Sec. III. In Sec. 



IV, it presents the formula in final integration to obtain 
the transverse momentum distribution of J/V' produc- 
tion. Sec. V. is devoted to the description about the 
calculation of J/ijj polarization. The color factor treat- 
ment for all the calculated processes are given in Sec. VI. 
In Sec. VII, treatment of T is given. The numerical re- 
sults are presented in Sec. VIII. Finally, The conclusion 
and discussion are given in Sec. IX. 



II. THE LO CROSS SECTION OF J/ip 
HADRONPRODUCTION 

The related Feynman diagrams which contribute to the 
LO amplitude of the partonic process g{pi) + g{p2) 
J /4'{Pz) + giPi) are shown in Fig.[Tl while the others can 
be obtained by permuting the places of gluons. 

In the nonrelativistic limit, we can use the NRQCD 
factorization formalism to obtain the partonic differential 
cross section in n = 4 — 2e dimension as 



dg^ _ b^al\Rs{Q)\''[s''{s-lf +P{i-lf +e{u-lf 
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by introducing three dimcnsionless kinematic variables: 

. ^ (Pl +P2f £ ^ [Pl - PJ.? ^ ^ {Pl - Pif 

Am^ ' 4to^ ' 4m^ ' 

where -Rs(O) is the radial wave function at the origin of 
J/^ and the reasonable approximation Mj/^ = 2mc is 
taken. 

The LO total cross section is obtained by convoluting 
the partonic cross section with the parton distribution 
function (PDF) Gg(a;,/i/) in the proton: 

= J dxidx2Gg{xi,^if)Gg{x2,fJ.f)a^ , (3) 

where Hf is the factorization scale. In the following a 
represents the corresponding partonic cross section. 

HI. THE NLO CROSS SECTION OF J/ip 
HADRONPRODUCTION 

The NLO contributions to the process can be written 
as a sum of two parts: one is the virtual correction which 
arises from loop diagrams, the other is the real correction 
caused by radiation of a real gluon, or a gluon splitting 
into a light quark-antiquark pair, or a light (anti)quark 
splitting into a light (anti) quark and a gluon. 
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FIG. 1: Leading order Feynman diagrams for g + g J/x}i + g. 
The other five diagrams can be obtained by permutation the 
places of gluons. 



A. Virtual corrections 

There are UV, IR and Coulomb singularities in the cal- 
culation of the virtual corrections. UV-divergences exist- 
ing in the self-energy and triangle diagrams are removed 
by the renormalization of the QCD gauge coupling con- 
stant, the charm quark mass, charm quark and gluon 
fields. Here we adopt renormalization scheme used in 
ref. fi^]. For the charm quark mass, charm quark and 
gluon fields, the renormalization constant Zm, Z2 and 
Z3 are determined in the on-mass-shell(OS) scheme while 
for the QCD gauge coupling constant, Zg is fixed in the 
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FIG. 2: One-loop diagrams for gg — > J/^g. Group (a) and 
(b) are counter-term diagrams of the quark-gluon vertex and 
corresponding loop diagrams, Group (c) are the quark self- 
energy diagrams and corresponding counter-term ones. More 
diagrams can be obtained by permutation of external gluons. 



modified- mimmal-subtraction(MS) scheme: 
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where 7e is Euler's constant, /3o = — Cyi — ^^F^f is 
the one-loop coefficient of the QCD beta function and 
is the number of active quark flavors. There are three 
massless light quarks w, d, s, so nj—i. In S'J7(3)c, color 
factors are given by Tp = 5; C'f — |, Ca — 3. And /i^ is 
the renormalization scale. 

After having fixed the renormalization scheme, there 
are 129 NLO diagrams in total, including counter-term 
diagrams. They are shown in Fig. [2l and divided into 8 
groups. Diagrams of group (e) that has a virtual gluon 
line connected with the quark pair lead to Coulomb sin- 
gularity, which can be isolated by introducing a small 
relative velocity v = |pc — Pel- The corresponding contri- 
bution is also of 0(as) and can be mapped into the cc 
wave function. 
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is nonzero. Otherwise, It is to do the integration directly 
with Feynman parametrization for two-point tensor case, 
and to write the Lorentz structure with independent ex- 
ternal momentums and apply the Passarino-Veltman re- 
duction again for other cases. In the calculation of scalar 
integral, we first try to decompose the scalar integral 
into several lower-point ones when it's Gram determi- 
nant is zero, if it fails, then to do the integration directly 
with Feynman parametrization just like the treatment to 
scalar integral with nonzero Gram determinants. Above 
procedure, including both reduction and integration, are 
done by FDC automatically. 

In our calculation, there are total 86 scalar integrals in 
total: 

• 65 of the total 86 integrals, can be found in Ref. Q 
after including the permutation of s, t and u. But 
the explicit results for the three Coulomb singular 
five-point scalar integrals is not available in Ref. 0] . 

• The remaining 21 integrals are not listed in Ref. 

12 of them can be reduced to combination of some 
lower-point scalar integrals and needn't to be inte- 
grated directly. 

• Another 6 of them can be expressed by the 
following two integrals, C(j)i,p3,mc,mc,mc) and 
D{pi,p4,p3 -I- p4, 0, mc, TOc, mc), through permuta- 
tion of s, t and u, where A, B, C, D, E are defined 
exactly the same as in Ref. They can be writ- 
ten into a linear combination of another two scalar 
integrals as: 

C(pi,p3,mc,mc,mc) 
= ^C'(-P3/2,-p3/2+Pi,0,mc,mc) 

-fic(p3/2,-P3/2+Pi,0,m„mc), (6) 

D{pi,p4,p3 +p4,0,mc,mc,mc) 

= ^D{p3/2,p3/2 ~ p2,-p3/2 - p4,mc,mc,mc,mc) 

+ ^D{-p3/2,p3/2 - P2, -P3/2 - P4, trie, TOc, TOc, TTlc). 

But in our calculation, they are calculated indepen- 
dently, and above relationship can be used to check 
all three scalar integrals. 

• The remaining 3 scalar integrals can be expressed 
by B{pi,mc,mc) through the permutation of s, t 
and u. 

More details about these 86 scalar integrals can be found 
at FDC homepage [50|. 

By adding all diagrams together, the virtual correc- 
tions to the differential cross section can be expressed 
as 



The Passarino-Veltman reduction [4lj is adopted in 
the tensor decomposition when it's Gram determinant 



dt 



cx 2Re(M^M*'*), 



(7) 
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FIG. 3: Feynman diagrams for first three real correction pro- 
cesses, (a) is for gg Jf-ip + gg and (b) is for gg —> J/ip + qq 
while (c) is for gq{q) J/ip + gq{q)- More diagrams can be 
obtained by all possible permutation of gluons. 



where is the amplitude at LO, and is the renor- 
malized amplitude at NLO. is UV and Coulomb fi- 
nite, but it still contains the IR divergences: 
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And the total cross section of virtual contribution could 
be written as: 



V 



(10) 



dx 1 dx2 Gg (a; 1 , /z/ ) Gg (x2 , M/ ) o" 



B. Real corrections 

The real corrections arise from four parton level sub- 
processes: 

9{Pl)+9iP2) J/i^iP3)+9{P4)+9iP5), (11) 

9{Pl)+9{P2) J/^KP3) + <l{Pi)+QiP5), (12) 

9iPi) + q{q){P2) J/V'(P3) + 5(P4) + g(g)(P5).(i3) 

9{Pi)+9{P2) J/VXp3) + c(p4) + c(p5), (14) 




FIG. 4: Feynman diagrams for real correction process gg 
J/tp + cc. More diagrams can be obtained by reversing the 
arrow of charm quark line and/or exchanging the places of 
gluons. 



above processes (except gg J/ip + cc) generates IR sin- 
gularities, which are either soft or collinear and can be 
conveniently isolated by slicing the phase space into dif- 
ferent regions. We use the two-cutoff phase space slicing 
method [42] to introduces two small cutoffs to decompose 
the phase space into three parts. 

Real gluon emission brings soft singularities. A small 
soft cutoff 5s is used to divide the phase space into two 
regions according to that the emitted gluon is soft or 
hard. Then another small cutoff 6c is used to decompose 
the hard region into collinear and noncoUinear regions. 
Then the cross section of real correction processes can be 
written as 



(15) 



The hard noncoUinear part o-^'^ is IR finite and can be 
numerically computed using standard Monte-Carlo in- 
tegration techniques. The subprocess gg — *■ J/ip + cc 
consists of only hard noncoUinear part. 



soft 



We have neglected the contribution from a real cor- 
rection subprocess qq —^ J/ipgg, which is IR finite and 
tiny (it only contributes about 0.002% at pt=3 GeV and 
0.05% at pt— 50 GeV to the differential cross section). 
And Feynman diagrams for above processes are shown 
in Fig. [3] and Fig. 21 The phase space integration of 



It is easy to find that soft singularities caused by emit- 
ting soft gluons from the charm quark-antiquark pair 
in the S-wave color singlet J/ip are canceled by each 
other. Therefore only the real gluon emission subpro- 
cess in Eq. pT|) . where there could be a soft gluon emit- 
ted from the external gluons, contains soft singularities. 
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Suppose P5 is the momentum of the emitted gluon. If we 
define the Mandelstam invariants as s^- = (pi +Pj)'^ and 
tij = {pi—pjY, the soft region is defined in term of the en- 
ergy of P5 in t\\epi-\-p2 rest frame by < -E5 < 5s^/sr2/2. 
In this region, soft approximation can be made and the 
matrix element squared can be factorized as 



i,i=l,2,4 



[Pi ■ P5){P] ■ P5) 



with 



and 



T^WMf,...,^,...,, T'^(j)Mf^...,^, 



(16) 
(17) 

(18) 



where M^ . j,^ is the color connected Born matrix ele- 
ment. 

Meanwhile, if we parametrize the emitted giuon's n- 
dimension momentum in the pi + p2 rest frame as 



Pb = -£'5(1, . . . , sin^i cos 6*2, cos 6*1) , 



(19) 



the three-body phase space in the soft limit can also be 
factorized as 



drslsoft = dr2 
with 
d^ 



Til 



1 



47r 



S12J r(l - 2e) 2(27r) 



dS , (20) 



1 

TT \S12 / JO 

X sin-2'^6'2 ^6*2 , 



dE^El-^' sini-2^6li dOi 
(21) 



as given in ref. [4^ . After analytical integration over the 
soft gluon phase space, the parton level cross section in 
soft region can be expressed 



a ~ (7 



with 



r(l-e) /47r^2 



277 r(l - 2e) V si2 
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A| = 9, Af = 3 
and 
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18 In ,5s, 
(23) 
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(24) 



2. hard collinear 

The hard collinear regions of the phase space are those 
where any invariant (s^ or tij) becomes smaller in mag- 
nitude than ScSi2. It is treated according to whether 
the singularities are from initial or final state emitting 
or splitting in the origin. Subprocess in Eq. (fT2|l con- 
tains final state collinear singularities, and subprocess in 
Eq. (|13p contains initial state collinear singularities while 
subprocess in Eq. pT|) contains both. 

a. final state collinear For subprocesses in Eq. (jlip 
and (|12p . the final state collinear region is defined by 
< S45 < ^12^ As a consequence of the factorization 
derivation |43l . |4J] , the squared matrix element factorizes 
into the product of a splitting kernel and the LO squared 
matrix element. 



\Mr\^ \coll 



S45 



B\2 



(25) 



where 4' denotes the parton which splits into parton 4 
and 5 collinear pair and (z, e) are the unregulated 
{z < 1) splitting functions in n = 4 — 2e dimensions 
related to the usual Altarelli-Parisi splitting kernels [i^ 
with z denoting the fraction of the momentum of parton 
4' carried by parton 4. For z < 1 the n— dimensional 
unregulated splitting functions are written as Pij{z, e) = 
Pi-j (z) + ePij (z) with 

1 + z^ 





= Cf 




= -C 
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= 0, 








P',ai^) 


= — z 



1-z ' 



1 - z 



+ z(l-z) 



w + (1 - 



(26) 



Meanwhile, the three-body phase space in the collinear 
limit can also be factorized as 1421: 



dP 



3 coll 



dr. 



(47r) 



167r2r(l -e) 



Azdsi^[si^z{l - z)]-' . (27) 



Hence after integrations of z and S45, the parton level 
cross section in hard final state collinear region can be 
expressed as 



as r(l-e) fA^^^ll 



27rr(l - 2e) V S12 
I ^ li^l ^ ^a-^gg ^ ^g-.qq \ 



(28) 
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where Ai and Aq are 
^9-99 ^ 3 (11/6 + 2 In (5^) 

= 3[67/18-7rV3-ln^(5;-ln,5, (ll/6 + 21n5;)] 
Ar"' = -nf/3 

Ag^m = n//3(ln,5c-5/3) , (29) 
for subprocesses in Eq. (fTT|) and ([T^ . and 

Sl2 s ^ 



6i 



Sl2 + S45 



1 



(30) 



Thus the total cross section for real correction processes 
in hard final state collinear region can be written as: 



dx 1 dx2 Gg (a; 1 , /X/ ) Gg (x2 , /i/ ) o"/ 



HC 



(31) 



b. initial state collinear For subprocess in Eq. (|13p . 
the hard initial state collinear region is defined by < 
—t25 < ^cSi2- However for subprocess in Eq. pT|) . the 
hard initial state collinear region is defined if any of the 
following conditions is satisfied < —tij < (5cSi2, with 
i = 1,2 and j = 3,4. For convenience, suppose that 2 
and 5 are the partons involved in the splitting 2 ^ 2' + 5 
while 2' denotes an internal gluon. Following the similar 
way as in the final state collinear case, the squared matrix 
element can be written as 



-zt 



25 



-P2'2{z,e)\M 



B\2 



(32) 



where z denotes the fraction of parton 2's momentum 
carried by parton 2' with parton 5 taking a fraction (1 — 
z). And the three-body phase space in the collinear limit 
can also be factorized as: 



dr 



3 I coll 



= dr2 



[Any 



167r2r(l - e) 
The i25 integration yields 

-cli25(-^25)"'" 





dzdt 
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-{5cSl2y 
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If we write the total cross section of LO as 
da^ — dxidx2Gg{xi)Gg{x2)da^ , 



(33) 



(34) 



(35) 



where Gg{xi) is the bare PDF. And using above results, 
the three-body cross section in the hard initial state 
collinear region can be written as [1^ 



da, 



HC 



Gg{xi)G2{y)dyda^{zsi2,ti3,ti4) 
as r(l — e) / 47r/i: 



27rr(l - 2e) V S12 



1 



— 5-'P2'2iz,e)dz{l- zY 



xS{yz — X2)da;ida;2 



(36) 



Notice that a factor of 1/z has been absorbed into the 
fiux factor for the two-body subprocess, and the delta 
function used here ensures that the fraction of hadron's 
momentum carried by 2' is X2 ■ And one more thing that 
need to be cared is, S12 here is related to the square 
of the overall hadronic squared center-of-mass energy S 
by S12 — xiyS, but in the LO process the relation is 
•S12 = X1X2S. From now on, we take the latter definition, 
so that the replacement S12 2/512/2^2 should be made. 
After the y integration we have 



da, 



HC 



= Gg{xi)G2ix2/z)da 



as r(l-e) /4^^2 



-i)<5r^2'2(z,6)^ 



27rr(l - 2e) V si2 
^ ^ dxidx2 . 



(37) 



When all possible two-to-three subprocesses are consid- 
ered, there will be several contributions, corresponding 
to a sum over all possible parton 2. It can he 2 — g fol- 
lowed hy g ^ gg or 2 = q{q) followed by q{q) — > 9(9) <?■ 
The collinear singularity must be factorized and absorbed 
into the redefinition of the PDF, which is in general called 
mass factorization [il] . Here we adopt a scale dependent 
PDF using the modified minimal subtraction (MS) con- 
vention given by [i^ . 



Gb{x,fif) = Gb{x) + — 



r(l - e) / 47r^2 



27rr(l - 2e) 



^ dz 

— Pbb'iz)Gb'ix/z) . 



(38) 



Use this definition to replace Gg{x2) in the LO expression 
psp and combine the result with the hard initial collinear 
contribution (j37p . then the resulting 0(a .) e xpression for 
the hard initial collinear contribution is 14211 



da. 



HC 



da 



as r(l - e) / 47r^^ 



27rr(l-2e) V S12 

X <j Gg{Xl,Hf)Gg{X2,^if) + 



Ar{g gg) 



(39) 



-At^ig^gg) 



Gg{xi,iif)Gg{x2,tJ.f) ^da;ida;2 



with 

Gc{x,fif) = 

c' 

and 



dy 



Gc'{x/y,fif)Pcc'{y), (40) 



p^J{y)^p.Ay)^^(sc^^^]-p:,{y). (41) 

The soft collinear factors Af^ result from the mis- 
match in the z integrations. They are given by Aq'' = 



7 



In(si2/Atp and (5 ^ gg) = 6 ln(5s + (33 - 2n/)/6. 
For subprocess in Eq. ([T^ . the light quark(antiquark) 
can come from either initial hadrons, while for subpro- 
cess in Eq. initial coUinear may happen to either of 
the initial gluons, thus the cross section of hard initial 
coUinear regions can be written as 

(^"^ ^<^add+ / Cr'j^^Gg{Xl,^lf)Gg{x2,^if)dXlAx2, 



with 



"add 



Us r(l — e) / Att^ 



,2 \ « 



27rr(l - 2e) V si2 
Gg{xl,flf)Gg{x2,^J,f) + {xi <-> X2] 



(42) 

(43) 

da;idx2, 
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and 



as r(l — e) / Anfi 
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Arig ^ 55) 



+ Kid ^99) 



(44) 



C. Cross section of all NLO contributions 



The cross section of real correction processes in hard 
noncoUinear regions could be written as 



HC 



{99 ^ J/i^ + 99)+ ^""^ i99 ^ J/^ + qq) 

q—u,d,s,c 



dxidx2Gg{xi, ^j.f)Gg{x2, fJ-f) 



J2 »"^{ga^ J/iP + ga) 



a.—u,d,s,u,d.s 



Gg(xi,/i/)G'Q(x2,/i/) + {xi ^ X2) 



dxidx2, 



(45) 



Thus the cross section of all real corrections becomes 

<^ = ^add+^ + /('^ +^/ ) 

xGg{xl,^J,f)Gg{x2,^J.f)dxldx2■ (46) 
And the total cross section of NLO QCD correction is 



NLO _ HC I HC I V 



xld 



G 



a 



(47) 



with 



G + (T + CT + (7^ + (Tj j 



X Gg (x 1 , ^/ ) Gg (a;2 , M/ ) da; 1 da;2 . 



(48) 



It is easy to find that there is no IR singularities in above 
expression, for 2A^ + = and 2J^ + Af + Af + 
_|_ 2A^y{9 ^ 99) — 0- The apparent logarithmic 
8s and 8c dependent terms also cancel after numerically 
integration over the phase space. 



IV. TRANSVERSE MOMENTUM 
DISTRIBUTION 



dptdy) is introduced. Thus we have 



dxi dx2dtGg (xi , /!/ )Gg {x2 , ) 



dG 

'dt 



.do- 



= I JdxidptdyGg{xi, fj.f)Gg{x2: fif) — , (49) 



and 



dG_ 

dpt 



JdxidyGg{xi,fif)Gg{x2,Hf)^, (50) 



with 



Pi =2^1^(1,0,0,1), P2 = a;2^(l, 0,0,-1), 



uit = ^/Mj,,+pf, p3 = (mt cosh 2/,pt,0,TOtsinhy), 



Xt 



J = 



Xi\ 



^XiX2Pt 

2xi — xtey ' 

_2t + Xtey 
^ 2 - xte-y ■ 



ml - M2 



2:2 



2t + xi XfC^y 
2xi — Xtey ' 



(51) 



To obtain the transverse momentum distribution of where ^/S is the center-of-mass energy of pp{p) at Teva- 



J/ip, a transformation for integration variable (da;2dt 



tron or LHC, 1714 is the invariant mass of all the final 
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state particles except J/V', and y and pt are the rapid- 
ity and transverse momentum of J/?/; in the laboratory 
frame respectively. 



V. POLARIZATION 

The polarization factor a is defined as: 



dar/dpt + 2daL/dpt 



(52) 



It represents the measurement of J/tp polarization as 
function of J/ip transverse momentum pt when calcu- 
lated at each point in pt distribution. To calculate a{pt), 
the polarization of J/tp must be explicitly retained in the 
calculation. The partonic differential cross section for a 
polarized J/^p could be expressed as: 



d^x , , 

-df^'"'^^^ 



aij Pt ■ e(A) pj ■ e*(A), (53) 



where A = ri,r2,L. e{Ti), t{T2), e{L) are the two trans- 
verse polarization vectors and the longitudinal polariza- 
tion one of J/V'i and the polarization of all the other 
particles are summed over in n-dimensions. It causes a 
more difficult tensor reduction path than that with all 
the polarizations being summed over in the calculation 
of virtual corrections. It is found that a and aij are fi- 
nite when the virtual corrections and real corrections are 
summed up. Therefore there is no difference in the dif- 
ferential cross section da\/dt whether the polarization of 
J/V' is summed over in 4 or n dimensions. Thus we can 
just treat the polarization vectors of J/ip in 4-dimension, 
and also the spin average factor goes back to 4-dimcnsion. 

To make a cross check, we carry out another calcula- 
tion. Namely, we calculate the differential cross section 
"'add cr^^ with the the polarizations of all particles 
being summed up analytically. The results are numeri- 
cally compared with that obtained without summing up 
the polarization of J/ip- Moreover, to check gauge invari- 
ance, in the expression we explicitly keep the gluon po- 
larization vector and then replace it by its 4-momentum 
in the final numerical calculation. Definitely the result 
must be zero and our results confirm it. To calculate 
, only numerical computation is carried out and we 
only sum over the physical polarizations of the gluons to 
avoid involving diagrams which contain external ghosts 
lines. 



VI. COLOR FACTOR 



processes, color factors are orthogonalized and normal- 
ized. There are three color factors in amplitude level for 
real correction process g + g ^ J li> + 9 + 9 

_}_^j.rjiC4 JiCi/jnC2 rpC4 rpC2 rpCi rpC^l f54) 

^/5 ^ J' V 7 

_}_'^j.[j^C-4rpCirpC2J-^C^ 2^C4 ^^£5 /^C2 1 

where Ci are the color indices of the external gluons. For 
9 + 9 ^ J /ip + q + q, there is one color factor 



75 

6%/5 



(55) 



where ci, C2 and C4, C5 are the color indices of the external 
gluons and quark, anti quarks respectively. And g + q^ 
J /ip + g + q has almost the same color factor as above. 
For g + g J/ip + c + c, there are three color factors 



1 



2V66 




52V15 

where ci , C2 and C4, C5 are the color indices of the external 
gluons and c quark, anti c quarks respectively. 



VII. TREATMENT OF T 

The production mechanism of T at Tevatron and LHC 
is much similar to that of J/ip except that, color octet 
states contribute much less in T production according 
to the experimental data and LO theoretical predictions. 
We can apply the results of above calculation to the case 
of T by doing the substitutions: 



'j/ip Mr 



Mj 

Uf = 3 ^ rt/=4 



(57) 



Note that charm quark is treated as light quark as an 
approximation. It is not coincide with the definition of 
CTEQ6M PDFs used in the calculation. The mass of 
heavy quark is not zero in the definition of CTEQ6M 
PDFs. This approximation can cause a small uncertainty. 



There is just one color factor dciC2C4 for the LO pro- 
cess in amplitude level with ci , C2 and C4 being the color 
indices of the three gluons in the process. And it is the 
same for the virtual correction process that just only one 



color factor dn 



appears in amplitude level. For other 



VIII. NUMERICAL RESULT 

In our numerical calculations, the CTEQ6L1 and 
CTEQ6M PDFs [i^, and the corresponding fitted value 
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for as{Mz) = 0.130 and as{Mz) = 0.118, are used for 
LO and NLO predictions respectively. At NLO, we use 
as in two-loop formula as 

^ 1 /3i lnln(A.VA^g^) 

where (3i = 34Cy3-4(Ci. + 5CA/3)Ti.n/ is two loop co- 
efficient of the QCD beta function. For the heavy quark 
mass and the wave function at the origin, rric = 1.5 GeV 
and |i?^(0)|2 = 0.810 GeV^ are used for J/V", and mt = 
4.75 GeV and \Rs{0)\^ = 0.479 GeV^ are used for T. To 
choose the renormalization scale fj,r and the factorization 
scale /i/ in the calculations is an important issue and it 
causes the uncertainties for the calculation. We choose 
jjL = IJLr = P^f = \J (2TOg)^ -I- vl as the default choice in 
the calculation with mg being rric and for J ji]) and 
T respectively. The center-of-mass energies are chosen 
as 1.98 TeV at Tevatron and 14 TeV at LHC. The two 
phase space cutoffs 5s and 8c are chosen as (5^ — 10~^ and 
8c = 8slh^ as default choice. To check the independence 
of the final results on the two cutoffs, different values of 
8s and 8c are used, where 8s can be as small as 8s = 10^^. 
And the invariance is observed within the error tolerance 
of less than one percent. 

It is known that the perturbative expansion cannot be 
applicable to the regions with small transverse momen- 
tum and large rapidity of J/4> or T. Therefore, Pt > 3 
GeV are used for all the calculations. For rapidity cut 
at Tevatron, we choose the same cut condition as the ex- 
periments at Tevatron [H, [3^: \y\ < 0.6 for J/ip and 
\y\ < 1.8 for T. To follow the same cut condition used 
in Ref. [s^, we choose \y\ < 3 for all calculation at LHC 
and another calculation of J/ip production at Tevatron. 
All the cut conditions are explicitly expressed for each 
result. 

The dependences of the total cross section at the renor- 
malization scale fir and factorization scale /i/ are pre- 
sented in Fig.[SJ Since the contribution from the subpro- 
cess gg J/4>cc is less than 10% of the total result at 
NLO in the whole reg ion of /x, it gives almost same plot as 
the Fig. 3 in ref. [3y] which does not included the contri- 
bution. The results show that the NLO QCD corrections 
boost the total cross section by a factor of about 2 at 
the default choice of the scales fi — fio ~ \/ {2mc)'^ + Pt- 
one can find that the scale dependence at NLO is not 
improved for J/tp- 

In Figs. [HI [3 and [51 the pt distribution of J/4> and T 
is shown. It is easy to see that the contribution of NLO 
correction becomes larger as pt increases, and in high 
Pt region, the NLO prediction is 2-3 order of magnitude 
larger than the LO one. As already known, the contribu- 
tion from subprocesses gg — > J/ ipcc or gg — > Tbb, which 
is also of 0{as), is large at high pt region. In order to 
compare with results in ref. [s^ and also to see how large 
is the contribution, result excluding this contribution is 
shown in the figures as NLO~. And we could see from 
the figures that the contribution from gg Tbb in T 



NLO LHC 

104 _ _ _ LO LHC 

: „„„. NLOTevotron 
LO Tevatron 

lyj/J<3 ond p, >3GeV 




0,5 C,6 0,7 0,80,9 ' 2 3 'I 5 6 7 

FIG. 5: Total cross section of J/tp production at Tevatron 
and LHC, as function of the renorm alization and factorization 
scale with fir = fJ-f = fJ- and no = a/ (2mc)^ + Pt- 
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FIG. 6: Transverse momentum distribution of J/ip produc- 
tion at LHC (upper curves) and Tevatron (lower curves). 
NLO~ denotes result excluding contribution from subprocess 
gg J/^cc. 



production is less than that from gg — s- J /ipcc in J /tp 
case. 

The Pt distribution of J /ijj and T polarization factor 
a is presented in Fig. [5] and Fig. [TOl We can see in the 
figures that a is always positive and becomes closer to 
1 as Pt increases at LO, and this figure means that the 
transverse polarization is more than the longitudinal one 
and even becomes dominant in high pt region. But there 
is dramatical change when the NLO QCD corrections are 
taken into account. For J/V', ot is always negative and 
becomes closer to -0.5 as pt increases, this new figure in- 
dicates that the longitudinal polarization is always more 



10 









Ml C Tevntrnn : 




NLC" Tevairon T 


r ' 


LO Tevatron 


r 


ly,/J<0.6 







l\ o' J/+ (GeV) 

FIG. 7; Transverse momentum distribution of J/tp produc- 
tion at Tevatron. NLO~ denotes result excluding contribution 
from subprocess gg J/tpcc. 
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FIG. 8: Transverse momentum distribution of T production 
at LHC (upper curves) and Tevatron (lower curves). \yr\ < 
1.8 and \yr\ < 3 is taken for Tevatron and LHC respectively. 
NLO~ denotes result excluding contribution from subprocess 
gg Tbb. 



than the transverse one and even becomes dominant in 
high pt region. Meanwhile the J/ip polarization in sub- 
process gg — > J l%pcc is near zero. By including contri- 
bution of this subprocess, the total result shown in the 
left diagram of Fig.[5]is closer to the experimental result. 
For T, a varies from positive to negative and becomes 
closer to -0.6 as pt increases, this indicates that the lon- 
gitudinal polarization becomes more and more, and even 
becomes dominant in high pt region. The T polarization 
of subprocess gg — s- Tbb is also near zero. But from the 
total result in Fig. [TO] we can see that this subprocess con- 



tributes less than the corresponding one in case of J/V"- 
Also, we find that the contribution from light quarks af- 
fects the Pt distribution of polarization less than 10%. 
When compare the figures for J/ij] with those for T, we 
can see that they are very similar to each other except 
that a is higher and even becomes positive in lower pt 
region for T. It could be understand by extending the 
curves for J /i}) to lower pt , because for a certain pt value 
in T production, it corresponds a lower pt in J / ij: pro- 
duction by just considering the energy scale. 

B y co mparing the exp erimental measurements for J/?/; 
[33 , l48l | and for T [ssl |49| at Tevatron with above re- 
sults, we could see that, although NLO corrections can 
boost the transverse momentum distribution of J/ 1/) and 
T very much, it is still an order of mangitude smaller 
than the experimental data. The color octet channels 
are still needed to explain the pt distribution. Thus the 
NLO prediction for the polarization of direct J /ip and T 
via color singlet channel could not be used to compare 
with experimental data. 

We can write contribution of each channel as 



^2 4 



CI 



where the overall e dependent factor 



1 



(59) 



(60) 



and the term 1/(1 — e)^ is from the gluon spin average 
factor l/(n — 2). When all the contributions are summed 
up, we have Y^C^ = Q and = 0. Thus Ce comes 

back to 1 and we have our result as Yl^o- Table. [H 
Cq is given. It should be careful that the Af^^l^g gg) 
term has been put into the gg — *■ J/ip + gg channel even 
if it contains a term proportional to the number of active 
flavors nj. 



i 


process 


C,UlO"nb) 




fraction 


1 


gg J/Tpg 


0.4061±0.0006 


0.2174 


0.1056 


2 


gg J/ipgg 


2.47±0.04 


1.32 


0.64 


3 


gg -> J/tpqq 


0.133±0.001 


0.071 


0.035 


4 


gq -> J/ipgq 


0.582±0.001 


0.312 


0.152 


5 


gg J/ipcc 


0.2583±0.0003 


0.1382 


0.0672 


E 


pp J/^ + X 


3.84±0.04 


2.06 


1.00 



TABLE I: lists of contributions from each channel to the NLO 
total cross section of J/ip hadronproduction at Tevatron in the 
region pt > 3 GeV and \yj/^\ < 3. We have set fj,r = jj,/ — fJ,o. 
Corresponding result for is 1.8682 x 10^ nb. 



IX. CONCLUSION AND DISCUSSION 

We have calculated the NLO QCD corrections to the 
J/ip and T hadronproduction at Tevatron and LHC. Di- 
mensional regularization is applied to deal with the UV 
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FIG. 9: Transverse momentum distribution of J/tp polarization at Tevatron and LHC. NLO denotes result excluding contri- 
bution from subprocess gg J/tpcc and the process itself is represented by the unlabeled dotted line. 
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FIG. 10: Transverse momentum distribution of T polarization at Tevatron and LHC. NLO denotes result excluding contri- 
bution from subprocess gg Tbb and the process itself is represented by the unlabeled dotted line. 



and IR singularities in the calculation, and the Coulomb 
singularity is isolated by a small relative velocity v be- 
tween the quark pair in the meson and absorbed into 
the bound state wave function. To deal with the soft 
and coUinear singularities in the real corrections, the two- 
cutoff phase space slicing method is used. By summing 
over all the contributions, a result which is UV, IR and 
Coulomb finite is obtained. 

Numerically, we obtain a K factor of total cross sec- 
tion (ratio of NLO to LO) of about 2 for J/^j. The 
transverse momentum distributions of J/V' and T are 
presented and they show that the NLO corrections in- 
crease the differential cross sections more as pt becomes 
larger and eventually can enhance it by 2 or 3 orders 
in magnitude at pt — 50 GeV. It confirms the calcula- 
tion by Campbell, Maltoni and Tramontano f36j. The 
real correction subprocesses gg — s- J/iJjcc and gg — s- Tbb 
are also calculated and the results are in agreement with 
those of Ref. [MIIJ. 



The NLO contributions to J/il; polarization is stud- 
ied and our results indicate that the J/4' polarization is 
dramatically changed from more transverse polarization 
at LO into more longitudinal polarization at NLO. All 
the results can be directly applied to 4' production by 
multiplying a factor (CJ* )/{Ot). The NLO contribu- 
tions to T polarization is also studied and presented for 
the first time. Our results indicates that at NLO, the 
polarization of T decreases gradually from near 0.2 to 
-0.6 as Pt increases from 3 GeV to 50 GeV. Namely, the 
Pt distribution of the polarization status behaves almost 
the same as that for J/tp except that the NLO result is 
also transverse polarization at small pt range. Since the 
fact that contribution via color-octet states is much less 
in T production than that in J/ijj case, our new result 
for T polarization plays an important role in understand- 
ing the experimental data. And even though our calcu- 
lation results in a more longitudinal polarization state 
than the recent experimental result for J/tp [3^ and T 
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[331 at Tevatron, it raises a hope to solve the large dis- 
crepancy between LO theoretical predication and exper- 
imental measurement on J/tjj and T polarization, and 
suggests that the next important step is to calculate the 
NLO corrections to hadronproduction of color octet state 
J/^'^^^ and T*^*). By re-fixing the color-octet matrix ele- 
ments, we will see what an involvement of the NLO QCD 
corrections can induce for the polarization of J/'0 and T. 
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